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Abstract 

QiQ ' This letter is a generalization of previous results on gravitational waves 

. (GWs) from f(R) theories of gravity. In some previous papers, particular 

f(R) theories have been linearized for the first time in the literature. Now, 
the process is further generalized, showing that every f(R) theory can be 

■ linearized producting a third massive mode of gravitational radiation. In 
^ ' this framework, previous results are particular cases of the more general 

problem that is discussed in this letter. The potential detectability of such 
massive GWs with LISA is also discussed with the auxilium of longitudinal 
response functions. 

Recently, the data analysis of interfcrometric GWs detectors has been started 
(for the current status of GWs interferometers seefTJISlOIllIllSJIHlEllE]) and the 
scientific community aims at a first direct detection of GWs in next years. 

Detectors for GWs will be important for a better knowledge of the Universe 
and either to confirm or ruling out the physical consistency of General Relativity 
or any other theory of gravitation [9l [T0 | [Til [121 [T3 l [T4] . 

An early, but quit interesting, work on GWs from other theories of gravity 
has been wriitcn by Wagoner 



1 



In this letter, the production of a hypotetical massive component of gravi- 
tational radiation which arises from a general f(R) theory of gravity is shown. 
The presence of the mass could also have important applications in cosmology 
because the fact that GWs can have mass could give a contribution to the dark 
matter of the Universe. 

The first and simplest f(R) theory of gravity was proposed by Starobinsky 
1151 . who dicussed the action 

(1) 

The production and the potential detection of GWs from this theory has been 
analysed in [15]. In [T7], it has been also shown that, from this particular 
linearized theory, it is possible to obtain an oscillating model of Universe. 
Another example is the action 

S = J d i x^j{R + R- 1 )+£ m . (2) 

This action has been analysed in a cosmological context in [18], while the 
GWs case has been analysed in [19]. This kind of theory could, in principle, be 
connected with the Dark Matter problem too [20]. Criticisms on f(R) theories 
of gravity arises from the fact that lots of such theories can be excluded by 
requirements of Cosmology and Solar System tests [H] . It is important to 
emphasize that the theory of eq. |J5J), differently from the theory of the action 
{1} is not in conflict with such constrains [24] [25]. In this context, even the 
theory 

j4~. rr^t„ r?i+ 



S= I d 4 xV^f R i+e + C m , (3) 
which has been discussed in [25] , is quite interesting. 

Equations |T]), ([2]) and ([3]) arc particular choices in respect to the well known 
canonical one of general relativity (the Einstein - Hilbert action [5TJ[52]) which 
is 

S = J d 4 Xv ^R + C m , (4) 

where R is the Ricci scalar curvature. 
Now, we will analyse the general case, i.e. 

S = J d 4 x^gf(R)+C m , (5) 

where f(R) is a generic high order theory of gravity. 

Of course, the cases which have been analysed in [16] and in [19] are partic- 
ular cases of the more general case that we are going to analyse now. 
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As we will interact with gravitational waves, i.e. the linearized theory in 
vacuum, C m = will be put and the pure curvature action 

S = J d A x^—gf{R) (6) 

will be considered. 

By varying the action ((6|) in respect to (see refs. [HI El [19] f° r a parallel 
computation) the field equations are obtained (note that in this paper we work 
with G = 1, c = 1 and h = 1): 

f'{R)R llv - l -f{R)g»u - f'{R)-,w + g^af(R) = (7) 

which are the modified Einstein field equations. f'{R) is the derivative of 
/ in respect to the Ricci scalar. Writing down, cxlplicitly, the Einstein tensor, 
cqs. (J7]) become 

G„„ = jjL-{l gia ,\f(R) - f'(R)R] + f(R) wv - 9»uUf\R)}. (8) 

Taking the trace of the field equations ([5]) one gets 

3D f'(R) + Rf'(R) - 2f(R) = 0, (9) 
and, with the identifications |23j 

*-/'(JJ) and % -> ^R)-Rf'(R) (10) 
a Klein - Gordon equation for the effective $ scalar field is obtained: 

M>=^. (11) 

To study gravitational waves, the linearized theory has to be analyzed, 
with a little perturbation of the background, which is assumed given by a a 
Minkowskian background plus $ = ^Oii-e- we are linearizing into a background 
with constant curvature [THl [23] ■ We also assume <&o to be a minimum for V: 

V ~ -a5<f> 2 => — ~ m 2 <5$, (12) 

and the constant m has mass dimension. 
Putting 

(13) 

$ = $0 + (5$. 

to first order in h^ u and 5$, calling R^ vpa , R^v and R the linearized quan- 
tity which correspond to R^pa , R^v and R, the linearized field equations are 
obtained QHQI1Q3]: 
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(14) 

□/ij = m 2 hf, 

where 

»' = •;■ < 15 > 

Rfivpo- and eqs. (fH)l arc invariants for gauge transformations [T51 [T71 [T5] 

(16) 

then 

^ = - + ?7m^/ ( 17 ) 
can be defined, and, considering the transformation for the parameter e M 

□e„ = V, (18) 
a gauge parallel to the Lorenz one of electromagnetic waves can be choosen: 

d M V = 0. (19) 
In this way, field equations read like 

□ V = (2°) 

Uh f = m 2 h f (21) 
Solutions of eqs. ([20]) and ([21]) are plan waves [l6l fT71 [T9] : 



h^v = A^ u (-f) e^p{ip a x a ) + c.c. (22) 
hf = a( ~f?) exp(iq a x a ) + c.c. (23) 



where 



fc Q = (cj,~p>) uj=p=[f\ 



(0J m ,^f) LU m = a/to 2 +p 2 



(24) 



In eqs. (|20p and (|22[) the equation and the solution for the standard waves 
of General Relativity [21] [22] have been obtained, while eqs. ([2T]) and (|23[) 
are respectively the equation and the solution for the massive mode (see also 

[moil us]). 
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The fact that the dispersion law for the modes of the massive field hf is not 
linear has to be emphatized. The velocity of every "ordinary" (i.e. which arises 
from General Relativity) mode /i M „ is the light speed c, but the dispersion law 
(the second of eq. ((Ml)) for the modes of hf is that of a massive field which 
can be discussed like a wave-packet fTHl El HH]- Also, the group- velocity of a 
wave-packet of hf centered in ~jf is 

vB = —, (25) 

UJ 

which is exactly the velocity of a massive particle with mass m and momen- 
tum J?. 

From the second of eqs. (fM|) and eq. (f2"5)l it is simple to obtain: 



\/iv 2 - m 2 

v G = . (26) 

Then, wanting a constant speed of the wave-packet, it has to be [TBI [T71 [T5] 



y/(l-vh)u. (27) 



Now, the analysis can remain in the Lorenz gauge with trasformations of 
the type \3e v = 0; this gauge gives a condition of transversality for the ordinary 
part of the field: k^A^ = 0, but does not give the transversality for the total 
field hy, v . From eq. (fT7|) it is 

hpv = h^v - + V^hf. (28) 

At this point, if being in the massless case [T51 H71 [TS] , it could been put 

□e' 1 = 

(29) 

0„e* = -§+/»/, 

which gives the total transversality of the field. But in the massive case this 
is impossible. In fact, applying the Dalembertian operator to the second of eqs. 
I|29p and using the field equations (|20[) and (f2"Tj) it results 

Ue^ = m 2 h f , (30) 



which is in contrast with the first of eqs. ()29|) . In the same way, it is possible 
to show that it does not exist any linear relation between the tensorial field h^ u 
and the massive field hf. Thus a gauge in wich h^u is purely spatial cannot be 
chosen (i.e. it cannot be put h^o = 0, see eq. (f2"5|)') . But the traceless condition 
to the field h^ u can be put : 



□e^ 1 = 



f) F » — - h. 



(31) 
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These equations imply 



0" V = 0. (32) 
To save the conditions d ll h^ v and h = transformations like 

□e^ 1 = 

(33) 

V = 

can be used and, taking ~]f in the z direction, a gauge in which only An, 
A22, and A12 = A21 arc different to zero can be chosen. The condition h = 
gives An = —A22- Now, putting these equations in cq. (f2"5]h it results 

V,(t, *) = - 2)e<f + A x (t - z) e W + ^(i - o G «)v (34) 

The term — z)e^ + A x (t — z)e\%) describes the two standard polar- 
izations of gravitational waves which arise from General Relativity, while the 
term hf(i—VGz)r\ ll , v is the massive field arising from the generic high order f(R) 
theory. In other words, the function f'(R) of the Ricci scalar generates a third 
massive polarization for gravitational waves which is not present in standard 
General Relativity. Note that the line element (fM|) has been obtained in both 
of references [TB] and [TO] starting by the actions (JTJ) and ^ respectively. Here 
we have shown that such a line element, i.e. the presence of a third massive 
polarization, is proper to every f(R) theory of gravity. 

The analysis of the two standard polarization is well known in the literature 
[21 [31 EH 122 ■ For a the pure polarization arising from the f(R) theory eq. 
can be rewritten as 

h^(t - v G z) = h f (t - v G z)r)^ v (35) 
and the corrispondent line clement is the conformally flat one 

ds 2 = [1 + h f (t - v G z)]{-dt 2 + dz 2 + dx 2 + dy 2 ). (36) 

In [19] it has been shown that in this kind of line element the effect of the 
mass is the generation of a longitudinal force (in addition to the transverse one) 
while in the limit m — > the longitudinal force vanishes. 

Now, before starting the analysis, it has to be discussed if there arc fenomcnog- 
ical limitations to the mass of the GW [THl [51] ■ A strong limitation arises from 
the fact that the GW needs a frequency which falls in the frequency-range for 
both of earth based and space based gravitational antennas, that is the inter- 
val 10- 4 Hz < f < WKHz [H H E M H H [H I2Z]. For a massive GW, from 
HUH EE EH it is: 

2tt/ = lo = Vm 2 +p 2 , (37) 
were p is the momentum. Thus, it needs 
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OeV < m < 1CT 1 W. (38) 

A stronger limitation is given by requirements of cosmology and Solar System 
tests on extended theories of gravity. In this case it is [23J 

OeV < m < 10~ 3S eV. (39) 

For these light scalars, their effect can be still discussed as a coherent GW. 

The frequency-dependent response function, for a massive mode of gravi- 
tational radiation, has been obtained in P]|] for the particular case f{R) = 
R + R~ l . Here, the computation will be performed with another treatment, 
showing that the longitudinal response function which has been found in [115] is 
proper to every f(R) theory of gravity, and the results will be applied to LISA, 
following the advice in [23J ■ 

Eq. ([3"o| can be rewritten as 

( | )M | )2 _ ( * )2 _ ( | )2 = _J_ m 

where r is the proper time of the test masses. 

From eqs. (|31)|) and (|30)) the geodesic equations of motion for test masses 
(i.e. the beam-splitter and the mirrors of the interferometer), can be obtained 







d 2 x 
dr 2 

& = o 

(41) 

d 2 t 1 dtjl+hf) 

Tr 1 2 (l + hf) 2 

d 2 z 1 d z (l+hf) 

dr 7 2 (l + hf)' 2 ' 

The first and the second of eqs. (|3Tj) can be immediately integrated obtaining 

Ci = const. (42) 



dx 



^ = C" 2 = const. (43) 
dr 

In this way eq. (j40]l becomes 

(-)»-(-)>= (44) 

If we assume that test masses are at rest initially we get C\ = Ci = 0. Thus 
we see that, even if the GW arrives at test masses, we do not have motion of 
test masses within the x — y plane in this gauge. We could understand this 
directly from eq. (|36|) because the absence of the x and of the y dependences in 
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the metric implies that test masses momentum in these directions (i.e. C\ and 
C2 respectively) is conserved. This results, for example, from the fact that in 
this case the x and y coordinates do not esplicitly enter in the Hamilton-Jacobi 
equation for a test mass in a gravitational field [2]. 

Now we will sec that, in presence of the GW, we have motion of test masses 
in the z direction which is the direction of the propagating wave. An analysis 
of eqs. (|4ip shows that, to simplify equations, we can introduce the retarded 
and advanced time coordinates (u, v): 

u = t — vqz 

(45) 

V = t + VqZ. 



From the third and the fourth of eqs. (|4Tj) we have 
d du d v [l + hf{u)\ 



drdr (l + h f (u)) 2 
This equation can be integrated obtaining 

du 



0. (46) 



dr * < 47 > 



where a is an integration constant. From eqs. (|44|) and (|47|) , we also get 

dv (3 



dr 1 + hf 

where /? = ^ , and 



(48) 



t = j3u + 7, (49) 

where the integration constant 7 correspondes simply to the retarded time 
coordinate translation u. Thus, without loss of generality, we can put it equal 
to zero. Now let us see what is the meaning of the other integration constant 
(3. We can write the equation for z from eqs. ([47|) and (j48]l : 



dz 1 , B 2 „. ,„„. 



When it is hf = (i.e. before the GW arrives at the test masses) eq. (|5t 
becomes 

I - > 2 - »■ ( 51 > 

But this is exactly the initial velocity of the test mass, thus we have to choose 
= 1 because we suppose that test masses are at rest initially. This also imply 
a = 1. 

To find the motion of a test mass in the z direction we see that from eq. 
(|49|) we have dr = du, while from eq. (j48]l we have dv = j^j- Because it is 
vqz = we obtain 



dz = ^TTh f - dul (52) 



which can be integrated as 

* = *a + 5fe/(i$7 -<*«) = 

z n --±- r*~ 1,GZ hf(u) du 



(53) 



2u G J- 

where zq is the initial position of the test mass. Now the displacement of 
the test mass in the z direction can be written as 



Az - Z - Zn - 3- ft-vezo-voAz h f (u) , 



^ 1_ rt-vczo h f (u) ^ 

2vg J— oo 

We can also rewrite the results in function of the time coordinate t: 

x{t) = x 

y(t) = yo 



(54) 



(55) 



T(t) = t - V GZ {t), 

Calling I and L + 1 the unperturbed positions of the beam-splitter and of the 
mirror and using the third of cqs. (|55[) the varying position of the beam-splitter 
and of the mirror are given by 

*Bs{t)=l-^tZ Gl ^d{u) 

(56) 

But we are interested in variations in the proper distance (time) of test 
masses, thus, in correspondence of eqs. (|56[) . using the fourth of eqs. (|55[) wc 



r M (t) = t - vqL -v G l -I tZ GiL+l) J^d(u). 
Then the total variation of the proper time is given by 



(57) 



1 ft-v G (L+l) U ( u ) 

A r(t) = r M (t) - T BS (t) = v G L - - / - - d(u). (58) 



2Jt-val l + h f (u) 
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In this way, recalling that in the used units the unperturbed proper distance 
(time) is T = L, the difference between the total variation of the proper time 
in presence and the total variation of the proper time in absence of the GW is 

1 ft-v a (L+l) h f (u) 

Sr(t) ee Ar(t) L = -L(v G + !)-■= " h i (59) 

This quantity can be computed in the frequency domain, defining the Fourier 
transform of hf as 

/oo 
dt hf(t)cxp(iut). (60) 
-oo 

and using the translation and derivation Fourier theorems, obtaining 
St(u) = L(l - v G ) exp[iwL(l + v G )] + 2ulL{ ^_ 1)2 

[exp[2iuL](v G + l) 3 (-2i + luL(v g - 1) + 2Lexp[iu>L(l + v G )] (61) 

{6iv G + 2iv% -uL + uoLv%) + L(v G + l) 3 (-2i + uL(v G + l))]h R . 
A "signal" can be also defined: 



5(c) ee = (1 - v%) exp[iu,L(l + vg)] + 2 ^ ( 4_ 1)a 

[ex.p[2iujL} (v G + l) 3 (-2i + tuL(v G - 1) + 2 exp[iwL(l + v G )} (62) 

(6iv G + 2iv G -loL + uoLv G ) + (v G + l) 3 (-2i + uL(v G + l))]h R . 
Then the function 

Ti(u) ee (1 - v 2 G ) exp[iwL(l + v G )] + 2ujL{ ^ g _ iY 2 

[exp [2iujL] (v G + lf(-2i + ljL{v g - 1) + 2 exp[iwL(l + v G )] (63) 
(6iv G + 2iv G ~ujL + ujLv g ) + (v G + 1) 3 (— 2i + ujL(v g + 1))], 

is the response function of an arm of the interferometer located in the z-axis, 
due to the longitudinal component of the massive gravitational wave arising from 
the high order gravity theory and propagating in the same direction of the axis. 
This longitudinal response function is proper to every f(R) theory of gravity 
and is quite interesting for f(R) theories of gravity which are not banned by 
requirements of Cosmology and Solar System tests [231 [55] . 

For vq — * 1 it is T;(a>) — > 0. Such a response function has been obtained in 
fH)] too, but with a different kind of analysis. 

In figures 1 and 2 are shown the response functions ([6"3")1 for an arm of LISA 
(L = 5 * 10 e Km) [5S[ [57] for v G = 0.1 (non-relativistic case) and v G = 0.9 
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Figure 1: the longitudinal response function (|63[) of an arm of LISA for vq = 0.1 
(non-relativistic case) 
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Figure 2: the longitudinal response function (|63p of an arm of LISA for vg = 0.9 
(relativistic case) 
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(relativistic case) . We see that in the non-relativistic case the signal is stronger 
as it could be expected (for nnOwe expectT/(w) — > 0). 

It is very important to emphasize that, differently from the response func- 
tions of massless gravitational waves, this longitudinal response function in- 
creases with frequency, .i.e , the presence of the mass prevents signal to drop off 
the regime in the high-frequency portion of the sensitivity band. Thus, consider- 
ing such a high-frequency portion of the sensitivity band becomes fundamental 
if LISA would detect massive GWs arising from f(R) theories of gravity which 
are not banned by requirements of Cosmology and Solar System tests [231 , 
like, for example, the two theories arising from the actions @ and ([3]). 

Conclusions 

This letter has been a generalization of previous results on gravitational waves 
(GWs) from f(R) theories of gravity. In some previous papers, particular f(R) 
theories have been linearized for the first time in the literature. Now, the process 
has been further generalized, showing that every f(R) theory can be linearized 
producting a third massive mode of gravitational radiation. In this framework, 
previous results are particular cases of the more general problem that has been 
discussed in this letter. The potential detectability of such massive GWs with 
LISA has been also discussed with the auxilium of longitudinal response func- 
tions. 
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